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ABSTRACT 

After the successful measurement of the mixing angle #13, the determination 

of the neutrino mass hierarchy has become a priority for future neutrino ex- 

f^) ' periments. We propose a conventional v^ beam with neutrino energies in the 

C^ • range 2-8 GeV aimed at a Mton underwater detector at the "magic" baseline 

?-h ' of 2600 km. In this constellation it is sufficient to distinguish (v^ induced) 

^"H track-like interactions from cascade-like interactions with moderate purity to 

^S ' determine the mass hierarchy. 

m i 

1. Introduction 

All three mixing angles and both mass square differences of the neutrino mass 
Oh 



Eigenstates are known after the measurement of the mixing angle #13 [B 12]. Global 
fits of all experimental input [3] provide a coherent picture of the oscillation param- 
eters. Among the yet unknown features of the oscillation scheme, the neutrino mass 
hierarchy (MH) is considered to be in reach for the next generation of experiments. 
j> \ It has been proposed to determine the MH by measuring atmospheric neutrinos [I] in 

planned low-energy extensions of existing neutrino telescopes [5j E]. The distinction 
of the two hierarchy hypotheses with atmospheric neutrinos is challenging due to can- 
cellation of contributions from neutrinos and anti-neutrinos and the effect is further 
attenuated by the finite energy and angular resolutions of the considered detectors. 
Both problems are avoided by counting beam related events of a specific flavour. 
Such a concept has been recently proposed [7J: muon counting from v^ interactions 
in very long baseline beams (L > 6000 km) should allow MH determination in Mton 
underwater/ice detectors. However, the construction of a steeply inclined beam-line 
has never been performed and is technically challenging and costly. Here, we propose 
instead to point a beam with a baseline of about 2600 km to a Mton underwater 
detector and to count v P interactions. 
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2. Oscillation Probabilities 

Oscillation probabilities are calculated in a full three flavour scheme using the 
Globes package [8] . For a given set of neutrino parameters for normal mass hierarchy 
(NH), a convention is needed to transform them into inverted hierarchy (IH). We 
follow the approach proposed in Ref. [3], defining an average "large" mass square 
difference AM 2 which controls atmospheric neutrino oscillations 

AM 2 = - (Am] 2)M + Aml hNH ) = - {Am 2 23JH + Am 2 13JH ) (1) 
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Figure 1: Neutrino oscillation probabilities P(v^ 



), P{y^ — > i^) and P{v^ — > f r ) for a baseline 



of 2600 km and oscillation parameters from a global fit as a function of neutrino energy. The group 
of red lines is for NH, blue for IH, <j>cp is varied in steps of 30° between 0° and 330°. The special 
values <ficp = 0°, 180° (CP conservation) are indicated in magenta, tficp — 90°, 270° (maximal CP 
violation) in orange on the left plot. 

with Am^ = m 2 — m%. IH is defined as the sign change of AM 2 for a given set of 
parameters in NH. This leads to the following transformation rules for the related 
mass square differences: 



Am?3,j/r 



Am 32,NH + Am 21 



^ m 31,NH 



Am 2 21 



(2) 
(3) 



Figure [T] shows the resulting oscillation probabilities Piy^ -> v e , u^, u T ) with the best 
fit parameters from Ref. [3] for a baseline of 2600 km. Most of the neutrino path will 
be in the outer Earth mantle with a maximal depth of 134 km. A constant density 
of 3.3 g/cm 3 is used for the calculation of the oscillation probabilities as given in the 
PREM model [9] for the outer mantle. 

For P{Pfj, — > u e ) a significant difference between both hierarchies is observed. 
The variation of the CP-phase <pcp (different lines of the same color) leads instead 
only to moderate changes of the oscillation probabilities. For E v = 3.5 GeV and IH 
P{V/i ~ * v e) is strictly independent from 4>cp- In the range 3 GeV< E u <8 GeV there 
is no overlap of the two CP-bands for the two hierarchies which allows to determine 
MH by counting v e events. This special feature has been noted in the past [10J and 
led to the label "magic" for baselines in the range 2500-2600 km. The large value of 
the mixing angle #13 ~ 9° leads to peak values of 13% for Piy^ — > v e ) in NH, allowing 
to detect a sizable sample of v e events in a suitable detector. 

P{ u n ~^ v y) depends only weakly on MH and not on <pcp- The first vacuum 
oscillation minimum at 5 GeV is the dominating feature in the shown energy range. 
Counting v^ events can serve as a flux normalisation or it could be used to improve 



the measurement of the atmospheric oscillation parameters AM 2 and #23 . 

3. Cross Sections 

The total cross sections for v^ and v^ charged current (CC) interactions are taken 
from Ref. [11] in the parton scaling approximation: 

a™{E v ) = 0.68-(£^GW)l(T 38 cin 2 (4) 

ag c \E u ) =0.3A-(E u /GeV)10~ 38 cm 2 . (5) 

Deviations from this linear behaviour due to quasi-elastic or resonant interactions are 
ignored. Their contribution would not alter the result of this study in a significant 
way. The relevant neutrino energies are significantly larger than m^ and m e therefore 
from flavour universality a^ c = of 7 and a^ = Cp . However for v T CC interactions 
the mass of the r-lepton cannot be neglected. We use the calculation from Ref. [12] . 
For a threshold energy E = 5 GeV and E < E v < 30 GeV we find the following 
simple parametrisation: 



c^) = 0.29bg(|j)a£ o (^ 



(6) 



The cross sections for neutral current (NC) interactions for all flavours are approxi- 
mated as: 

rt) =1< C (E») (7) 

<$°{E V ) =\o°?{E v ). (8) 

Whereas the sum of the oscillation probabilities J2 a P(^ ~^ a ) equals unity due to 
the unitarity of the flavour mixing matrix, the cross section weighted sum 

(normalised at Eq — 1 GeV) may help identifying the optimal energy range to separate 
the two MH hypotheses. P° is shown on the left plot of Figure[2]and can be interpreted 
as the event rate per neutrino energy seen in a detector with an energy independent 
detection efficiency for a pure z/ M flux flat in energy. No distinction is made neither 
between NC and CC interactions nor between different flavours. Nevertheless a clear 
separation of the two MH hypotheses is observed for energies above 3 GeV. It can be 
attributed to the kinematical suppression of v T CC interactions. The size of the effect 
is quantified in the right plot of Figure [2] which shows the ratio between IH and NH 
integrals from the left plot (for different values of <ficp) taken between an assumed 
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Figure 2: Left: Summed oscillation probabilities P(fi — > a) weighted by cross sections (red : NH, 
blue : IH) Right: Ratio of Integrals Nih/Nnh from left plot between E min =2.5 GeV and E max . 

threshold energy of 2.5 GeV and a variable maximal energy E max . An optimal value of 
E m ax = 6 GeV can be read from the figure, which yields a suppression of the IH event 
rate of 11-14% compared to the expected rate for NH. Extending the energy range 
to higher values reduces the relative size of the separation of the MH hypotheses. 
It can be concluded that a measurement of the mass hierarchy will be possible even 
without any flavour tagging capabilities using a neutrino beam in a limited energy 
range of 2-6 GeV and a large detector which can reliably count beam related neutrino 
interactions. Nevertheless flavour tagging methods are discussed below and they are 
used to improve the significance of the measurement. 



4. Neutrino beam 

The Institute of High Energy Physics (IHEP) in Protvino near Moscow (located at 
54°52'iV, Z7°11'E) hosts the U70 proton accelerator [T3] which provides protons with 
energies up to 70 GeV. It is operational since 1967 and had been the world largest 
proton accelerator at its time of commissioning. The beam provides 1.7- 10 13 protons 
per pulse at a repetition rate of 0.1 Hz. A planned upgrade to a total beam power 
of 0.45 MW will yield 2.2 ■ 10 14 protons per pulse at a repetition rate of 0.2 Hz [13J. 
This will allow for the accumulation of N pot = 10 21 protons on target after three 
years of operation assuming a running time of 6 months per year at a moderate beam 
efficiency of 50%. In the following, event numbers are calculated for N pot = 10 21 . 

In the past, a neutrino beam was provided to several experiments, among them 
the bubble chamber SKAT p3]. Secondary hadrons were produced in an Aluminum 
target and focused by parabolic lenses. Neutrinos were produced in a 140 m long 
decay line. The bubble chamber was situated 270 m behind the target, and at a 
distance of Iskat = 245 m downstream the beginning of the decay tube. This value 
will be used in the following to scale the beam intensity to the remote location. The 
neutrino fluxes d$ v /dE Vl as they were delivered to the SKAT experiment for focusing 
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Figure 3: Neutrino Flux d$ v /dE v per proton as seen by the SKAT bubble chamber. The different 
lines correspond to the different neutrino flavours. 

of positively charged hadrons and for all four flavours present in the beam, are shown 
in Figure |3j A parametrisation of the beam intensity as function of the neutrino 
energy at the SKAT detector has been obtained from Ref. [HI E] and will be used 
in the following. As seen from Figure [3j a very clean v^ beam was provided with 
energies dominantly in the range 2-8 GeV. Contaminations from other flavours were 
on the sub-percent level. 

To provide a beam to a detector in the Mediterranean Sea at a distance of 
Ilbl ~ 2600 km a new beam-line is needed. It would point in the southwest di- 
rection from the proton accelerator ring, an area which is currently not obstructed 
by buildings. A moderate downward inclination of 11.7° is needed. A 300 m long 
decay tunnel could be envisaged which yields 50% more neutrinos compared to 140 m 
used in the past. It would end 60 m below ground. A possible near detector might 
be located 500 m behind the target at a depth of 100 m. The calculation of event 
numbers below will be based on the traditional beam [H] scaled up by a factor 1.5 
to account for a decay tunnel of 300 m. 

5. Detector 



The ANTARES detector [16] is a deep sea neutrino telescope, which operates 
successfully in the Mediterranean Sea at (42°48'A, 6°10'E). A neutrino beam from 
Protvino to this location would result in a baseline of 2588 km. Recently, the 
ANTARES Collaboration published a measurement of atmospheric neutrino oscil- 
lations [T7], demonstrating the capability of the device to detect and measure neutri- 
nos with energies as low as 20 GeV. A multi-cubic-kilometer detector KM3NeT [18] 
is planned as a future neutrino telescope in the Mediterranean Sea. One possible site 



would be close to the existing ANTARES detector. An alternative site in the Ionian 
Sea off the Sicilian coast results incidentally in an identical baseline (within 1%) for 
a beam from Protvino. Therefore the calculation, presented here, holds for both site 
options. 
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Figure 4: Left: Effective mass M e ff for the detection of v^ CC interactions as function of E v for 
events with the interaction vertex inside the instrumented volume. A successful track reconstruction 
including a condition on the likelihood of the track fit is required. Right: Flavour misidentification 
probability as function of E v . The same function is assumed for the probability to misidentify tracks 
as cascades and vice versa. 



Currently the KM3NeT Collaboration carries out a feasibility study to determine 
the physics reach for MH determination of a densely equipped detector based on 
about 20% of the final budget [6J. About 1000 optical modules would be used to 
instrument a water mass of around 2 Mtons. The detector is supposed to be sensitive 
to neutrinos with energies down to few GeV. Figure H] shows a preliminary result of 
this study [6]: the effective mass M e ff for v^ charged current interactions occurring 
inside the instrumented volume is given as function of the neutrino energy. Here a 
trigger condition, a likelihood based track reconstruction and a quality criteria for 
the goodness of the track fit are applied. 

It is assumed that the considered detector is able to distinguish "track" and "cas- 
cade" event signatures. A "track" requires the presence of a long muon trajectory. As 
the muon range increases linearly with energy whereas the longitudinal extension of 
cascades remains essentially unchanged for the considered energies, the distinguisha- 
bility of these two topologies will be a function of the neutrino energy. For the event 
rate calculation below, we introduce e(E„) as the probability to misidentify a cascade 
as track, while the misidentification of a track as cascade will be called r](E v ). Both 
functions are parametrised by 



e(E v ) = r){E v ) = l/{E„/GeV); E v > 2GeV 



(10) 



which is illustrated in the right plot of Figure HI For E v = 2 GeV, e = rj = 0.5, 
which means the two topologies cannot be distinguished and the attribution of an 



event to one of them is random. A 5 GeV neutrino produces a muon with an average 
range of 15 m in Sea water which exceeds already by far the typical longitudinal 
size of a hadronic or electromagnetic shower. Correspondingly the misidentification 
probability is assumed to drop to 20%. For E v = 10 GeV e and r\ further decrease to 
10%. However, these low values of the misidentification probability are not exploited 
in this analysis due to the energy profile of the beam, which suppresses contributions 
of neutrinos with E v > 10 GeV. 



6. Signal Event Rates 

The rate of detected CC events of flavour a can now be calculated: 



dNa =N ( Iskat Y M eff {E v 
dE v pot \ l LBL ) m p 



a 



cc 



r act T . t- r pj 



dE tJ 



dE„ 



i ^cc 



d<S>p 
~d~E, 



d<S>u 
.p I ^p_ 



dE tJ 



;n) 



with m p the proton mass and the abbreviation Pp a = P{y$ — ¥ u a ) and Pj^ = P{vp — > 
u a ) for the oscillation probabilities. All four initial-state neutrino flavours which are 
present in the neutrino beam are taken into account. The major contribution for all 
final-state flavours come from the dominant i^-beam component. 

Information on reconstruction efficiencies for v e and v T CC events are not available 
yet. Contained v e and v^ CC events will on average deposit the same amount of 
energy in the detector, therefore the same neutrino energy dependent effective mass 
M e ff(E u ) from Figure H]is assumed for both. The same efficiency function is also used 
for v T CC events, despite the fact that they release less energy in the detector due to 
the escaping neutrino(s) from the tau decay. As v T interactions are a background in 
the present analysis, this is a conservative approximation. 

Similarly the rate of detected NC events can be calculated 
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No dependence on oscillation parameters enters here. As on average 50% of the 
neutrino energy is transferred to the outgoing neutrino, the effective mass M e ff is 
evaluated at E v /2. 

The "track-signal" events are then simply given by dN^ ck /dE u = dN^/dE u , 
whereas the "cascade-signal" will be defined as dN c s ^ c j dE v = dN e /dE u . v T CC and 
NC interactions are considered as background. Figure [5] shows the expected signal 
event rates. The corresponding integral event numbers are summarized in Table [TJ At 
this stage we do not consider the misidentification between the two event topologies, 





Figure 5: Observed events for N pot = 10 21 as function of the neutrino energy. Left: Track- like events 
from u^ CC interactions. Right: Cascade-Like events from v e CC interactions. The solid blue 
histogram on the left plot is obtained by ignoring neutrino oscillations. Solid magenta histograms 
are for NH and various values of 4>cp, dashed lines indicate IH. 



i.e. e — r\ — 0. The v^ CC rate is suppressed by 60% due to neutrino oscillations 
compared to the no-oscillation hypothesis. Nonetheless a comfortable event sample 
of 10000 events can be detected for N pot = 10 21 with a statistical uncertainty of 1% 
whereas the two mass hierarchy hypotheses modify the expected event rate by 3%. 



Channel 


Tracks NH 


Tracks IH 


Cascades NH 


Cascades IH 


No oscil 


26315 


— 


Signal 


10317 


10015 


1366-1876 


397-597 


NC 


227-231 



245-248 



1076-1098 
4732 


1163-1176 

4732 


BG Total 


227-231 


245-248 


5807-5830 


5895-5908 


Total 


10543-10548 


10260-10263 


7196-7683 


6304-6492 



Table 1: Event numbers for N pot = 10 21 in the track and cascade channel for both mass hierarchy 
schemes and varying 4>cp values with a perfect separation of track and cascade signatures. 

The right plot of Figure |5] shows the event rate of v e CC events. For NH 1621±255 
events are expected, where the uncertainty is due to the unknown CP-phase. This 
has to be compared to 497 ± 100 events for IH. The statistical separation of both 
samples is better than 20<r. The largest effect is seen for neutrino energies from 3- 
8 GeV, as expected. Despite the fact that Figure [5] illustrates the E v distribution of 
the selected event samples, no assumption about the energy determination is needed 
here to determine MH. However, a moderate energy resolution will certainly increase 
the significance of the MH hypothesis test. 

7. Background and Purity of Event Selection 

The measurement will be complicated by background. Contributions from at- 



8 



mospheric neutrinos and misreconstructed down-going atmospheric muons can be 
ignored. A pulsed beam with a typical duty cycle lower than 10 -6 allows to safely 
discard these events. The beam itself is a source of background events. v T CC events 
with a muonic r-decay produce "track" events, whereas NC events and v T CC events 
with a non-muonic r-decay have a genuine cascade signature. These two contributions 
are added in Table [1] and they are the only background contributions, if we assume 
a perfect flavour tagging mechanism to separate track-like from cascade-like events. 
Whereas the sample of track-like events is only marginally affected by the addition of 
v T events, the cascade event sample is now dominated by NC events which contribute 
about three times as much as the v e signal events. 

In a real detector it will not be possible to separate track-like and cascade-like 
events with 100% efficiency as discussed in Section 4. The probability to misidentify 
a cascade as track e(E u ) and the misidentification probability of a track as cascade 
T)(E U ) as introduced in Section 5 will be used now. With these two quantities the 
total background for the two event samples can be written 



J ]\Ttrack J ]\Tcasc 



dE v 
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(13) 



BR TfJL stands for the muonic branching ratio of the tau decay (17.4%). The total 
number of observed events in each channel is given by adding the (reduced) number 
of signal events and the background contribution from Equation | 



dNiT k 
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rf AT case 
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(14) 
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The resulting rates are shown in Figure |6] and quoted in Table [2j The signal contribu- 



Channel 


Tracks NH 


Tracks IH 


Cascades NH 


Cascades IH 


No oscil 


26315 


— 


Signal 


8990 


8735 


1134-1547 


350-519 


Misreco 

v T 

NC 


232-329 

324-332 

1092 


47-79 

351-355 

1092 


1326 

978-998 

3640 


1280 

1057-1068 

3640 


BG Total 


1655-1745 


1494-1522 


5944-5964 


5977-5988 


Total 


10645-10736 


10229-10257 


7099-7491 


6338-6496 



Table 2: Event numbers for N po t = 10 21 in the track and cascade channel for both mass hierarchy 
schemes and varying <ficp values. 



tions are reduced by 15-20%. The track-like sample is still dominated by signal. The 
different background channels add up to 15-18% of the overall rate. The situation for 
the cascade-like events does not change very much compared to the numbers given in 
Table [TJ Some backgrounds are reduced due to migration into the track-like channel 
(e.g. NC), others are increased. The total number of background events is almost 
unaffected by the choice of MH (and by 4>cp) as the MH dependence of different 
backgrounds has the tendency to cancel. This leaves the event difference of the signal 
part almost unaffected. The event rates between the two MH hypotheses differ now 
by 9-18% with a statistical uncertainty of 1.2%. The statistical significance of the 
MH hypothesis test is still better than la and it remains at the level of 3<r even 
when adding an additional systematic uncertainty of 3-4% (depending on the true 
value of the CP phase) for the determination of the total cascade event rate. The 
knowledge of the detector performance, water parameters, neutrino cross sections, 
oscillation parameters and the monitoring of the neutrino flux contribute to the sys- 
tematic uncertainty. The neutrino flux normalisation can be controlled by performing 
a complementary measurement of v^ CC events. Uncertainties of the oscillation pa- 
rameters will have been reduced by ongoing experiments and neutrino cross sections 
will have been measured with high precision by ongoing and planned short-baseline 
experiments by the time the here proposed experiment is running. The water in the 
abyss of the Mediterranean Sea is extensively studied in the ANTARES experiment. 
All these measurements will help reducing the systematic uncertainty. 
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Figure 6: Observed events for N po t — 10 21 as function of the neutrino energy. Left: Track-like 
events. Right: Cascade-Like events. In green the background from NC events is given. The blue 
histogram shows the total background including misreconstructed CC events. Magenta histograms 
are for the total event rate (signal plus background) for NH (solid) and IH (dashed) and various 
values of 0c p- 

Figure [6] illustrates the total event numbers as function of neutrino energy, detail- 
ing the different contributions. NC events are added to the E v /2 bins according to 
their lower light yield in the detector. Their main contribution is found for energies 
below 6 GeV. Background from CC events is instead mainly seen above 5 GeV. The 
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separation of the MH hypotheses is most pronounced in the range 4-8 GeV. 

8. Conclusion 

A neutrino beam from IHEP Protvino to a Mton detector installed in the abyss 
of the Mediterranean Sea has a baseline close to the "magic" value of 2600 km. A 
low energy "phase- 1" part of the future KM3NeT neutrino telescope could serve as 
target for such a beam. Counting of cascade-like events would allow a measurement 
of the neutrino mass hierarchy with a significance of 3<r for 10 21 protons on target 
and a systematic uncertainty of 3-4% for the event rate determination. Higher values 
of the significance can be reached by reducing the systematic uncertainty, considering 
a moderate capability to measure the neutrino energy of the signal events or an im- 
proved flavour tagging method. The proposed measurement would be complementary 
to an analysis of atmospheric neutrinos in the same detector. The combination of 
both measurements will possibly allow an unambiguous determination of the neutrino 
mass hierarchy. 
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